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New method for a continuous determination of the spin tune in storage rings 

and implications for precision experiments 
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A new method to determine the spin tune is described and tested. In an ideal planar magnetic 
ring, the spin tune - defined as the number of spin precessions per turn - is given by i/s = 7 • G (7 
is the Lorentz factor, G the gyromagnetic anomaly). At 970MeV/c, the deuteron spins coherently 
precess at a frequency of « 120 kHz in the Cooler Synchrotron COSY. The spin tune is deduced from 
the up-down asymmetry of deuteron-carbon scattering. In a time interval of 2.6 s, the spin tune was 
determined with a precision of the order 10“®, and to 1 x 10~®® for a continuous 100 s accelerator 
cycle. This renders the presented method a new precision tool for accelerator physics; controlling 
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the spin motion of particles to high precision is mandatory, in particular, for the measurement of 
electric dipole moments of charged particles in a storage ring. 

PACS numbers: 13.40.Em, ll.30.Er, 29.20.D, 29.20.dg 


The matter-antimatter asymmetry that emerges from 
the Standard Model (SM) of particle physics falls short 
by many orders of magnitude compared to the observed 
value [l|. Physics beyond the SM is thus required and is 
sought at high energies and by high-precision measure¬ 
ments at lower energies, for instance in the search for CP- 
violating electric dipole moments (EDMs). A non-zero 
EDM measurement would be a telltale sign of physics 
beyond the SM [^. In addition, EDM measurements of 
various systems would point towards the underlying ex¬ 
tension of the SM iil. 

Up to now, upper limits of hadronic EDMs have been 
determined for the neutron and the proton, but the 
latter only indirectly from a measurement on ^®®Hg Q. 
EDMs of charged hadrons are proposed to be measured 
in storage rings with a precision of 10“^® e-cm by observ¬ 
ing the influence of the EDM on the spin motion 
The high level of sensitivity is maintained only when the 
particle spins in the machine precess coherently for long 
periods of time (« 1000 s). In a series of recent investiga¬ 
tions at COSY 11, 1^, we studied how the spin coher¬ 
ence time of an ensemble of particles can be increased to 
many hundreds of seconds through sextupole corrections, 
bunching, and phase-space cooling of the beam 13 M- 

Another limiting factor of the storage ring approach 
to EDM searches, however, is controlling the spin mo¬ 
tion in the presence of small fluctuations of electric and 
magnetic fields in order to unambiguously determine the 
EDM signal. Consequently, the measurement described 
in this paper constitutes one cornerstone of storage ring 
EDM searches; viz. the first precise measurement of the 
spin tune during a complete accelerator cycle. 

The spin motion of a particle in the electric and mag¬ 
netic fields of a machine is governed by the Thomas-BMT 
equation 17-1^, extended to include the EDM 211, 


ds 

dt 


= s X ID 


+ ^edm) • 


( 1 ) 


Here, s'denotes the spin vector in the particle rest frame 
in units of h, t the time in the laboratory system, and 
Hmdm and Dedm the angular frequencies due to mag¬ 
netic dipole (MDM) and electric dipole moments (EDM). 
In the following, spin rotations due to EDMs, being 
many orders of magnitude smaller than those produced 
by MDMs, are neglected. It is convenient to define the 
spin motion relative to the momentum direction (^ . 
rotating with the angular velocity |Hcyc| = qB/imj). As 
a result, the spins of particles that orbit in an ideal pla¬ 
nar machine precess about the vertical magnetic field B 
relative to the momentum vector with the angular fre¬ 
quency Omdm = qGB/m, where q and m denote particle 


charge and mass, G is the gyromagnetic anomaly and B 
the magnetic field at a given point of the particle tra¬ 
jectory. Dividing |Hmdm| by the cyclotron angular fre¬ 
quency |Hcyc| yields the number of spin revolutions per 
turn, called the spin tune Vs 2^ 241. For a particle on 
the closed orbit in an ideal magnetic ring, the spin tune 
is thus given by 


lZs=jG. 


( 2 ) 


In a real machine, field imperfections, magnet misalign¬ 
ments, and the finite emittance of the beam lead to spin 
rotations around non-vertical axes and the spin tune de¬ 
viates from the one given in Eq. The most prominent 
one is the so-called pitch correction (25, 26|. 

The experiment was performed at COSY. A polar¬ 
ized deuteron beam of « 10® particles was accumulated, 
accelerated to the final momentum of 970MeV/c, and 
electron-cooled to reduce the equilibrium beam emit¬ 
tance. The beam polarization, perpendicular to the 
ring plane, was alternated from cycle to cycle using two 
vector-polarized states, = 0.57 ± 0.01 and p^~ = 
—0.49 ± 0.01, and an unpolarized state. The tensor po¬ 
larization pjj of the beam was smaller than 0.02. An rf 
cavity was used to bunch the beam during the ~ 140 s 
long cycle. After the beam was prepared, the electron 
cooler was turned off for the remaining measurement pe¬ 
riod of 100 s. 

An rf solenoid-induced spin resonance was employed 
to rotate the spin by 90° from the initial vertical di¬ 
rection into the transverse horizontal direction. Subse¬ 
quently, the beam was slowly extracted onto an inter¬ 
nal carbon target using a white noise electric field ap¬ 
plied to a stripline unit. Scattered deuterons were de¬ 
tected in scintillation detectors, consisting of rings and 
bars around the beam pipe (2^ , and their energy deposit 
was measured by stopping them in the outer scintilla¬ 
tor rings. The event arrival times with respect to the 
beginning of each cycle and the frequency of the COSY 
rf cavity were recorded in one long-range time-to-digital 
converter (TDC), i.e., the same reference clock was used 
for all signals. The number of orbit revolutions could thus 
be unambiguously assigned to each recorded event HI- 

In the following, we use a right-handed coordinate sys¬ 
tem, where the z-axis points in the beam direction, y 
upwards, and x sideways. The differential cross section, 
for scattering of purely vector-polarized deuterons with 
a vertical polarization component Pu = 0 off an unpolar¬ 
ized target, can be written as |28l l29| 


= ao{-d) 


1 - sincj) 


( 3 ) 
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Here, i7q(i9) denotes the differential cross section for un¬ 
polarized beam, i9 the polar scattering angle, 4> the az¬ 
imuthal scattering angle, and Ay{'d) the deuteron vector 
analyzing power. According to Eq. ([1]), Px(t) in Eq. ([3]) 
oscillates as 


Px {t) = Pi sin(Hst -K <f) , (4) 

where Hg = 2'Kf-ce^Va denotes the angular frequency of 
the horizontal spin precession, (p the phase, and = 
^Jp^ -f p% the magnitude of the in-plane vector polariza¬ 
tion. Because of the COSY straight sections, /rev differs 
from the cyclotron frequency /eye = Heye/(27r). 

In order to determine the spin tune from Eqs. ([3]) and 
a, and to cancel possible acceptance and flux varia¬ 
tions during the measurement, asymmetries are formed 
using the counts of the Up (U) and Down (D) detector 
quadrants. The quadrants are centered at (pu Ki 90° and 
(/d ~ 270°, covering polar angles from d = 9° to 13°, 
and an azimuthal range of A/)u ~ A(/d ~ 90°. An ex¬ 
pression for the event rate Rx of a detector quadrant 
X = (U or D) is obtained by integration over the solid 
angle, yielding 

Rx = Idt [ ax(d, (/))(T(d,(/)dH 
Jx 

= Idtopx (^1 - ^Pxit)A^^'^ . (5) 

Here, ax{'&,4’) denotes the combined detector efficiency 
and acceptance, I [s“^] the beam intensity, dt [cm“^] the 
target density, opx the integrated spin-independent cross 
section, and « 0.4 the weighted average analyzing 

power of the respective quadrants. 

It is not possible to determine the spin tune Vg from the 
observed event rates by a simple fit with Vg as a parameter 
using Eq. (0, because, at a detector rate of w 5000 s ^ 
and a spin frequency of /g = jr'gl • /rev ~ 0.16 • 750 kHz = 
120 kHz, only about one event is detected per 24 spin rev¬ 
olutions. Therefore, as described below, an algorithm is 
applied that maps all events into one oscillation period. 
It generates an asymmetry, largely independent of varia¬ 
tions of acceptance, flux, and polarization that oscillates 
around zero. For each event, the integer turn number n 
is calculated, using the event time compared to the time 
of the COSY rf cavity. Based on the turn number, the 
100 s measurement interval is split into 72 turn intervals 
of width An = 10® turns (each turn lasting « 1.3 ps). For 
all events, the spin phase advance ps = 27r|i/g |n is cal¬ 
culated under the assumption of a certain spin tune n®. 
In this analysis, we use the absolute value \vs\ because 
the asymmetry measured at one location in the ring is 
insensitive to the sign as well as to integer offsets of the 
spin tune. Each of the turn intervals is analyzed inde¬ 
pendently, and the events are mapped into a Itt interval, 
which yields the event counts N\j{pa) and Nd{ps) shown 
in Fig.[T](a). In order to obtain from Au(</?s) and Nd{ps) 



FIG. 1. (a): Counts N\j and Ad after mapping the events 
recorded during a turn interval of An = 10® turns into a spin 
phase advance interval of Itt. (b): Count sums Npp and 

differences of Eq. ([6| with € [0, 27r) using the 

counts N\j{pb) and Ad(v5s), shown in panel (a). The vertical 
error bars show the statistical uncertainties, the horizontal 
bars indicate the bin width. 


a sinusoidal wave form that oscillates around zero, four 
new event counts for the two quadrants {X = Uor D) are 
defined, 




Nx{Ps) ± Nx {Ps + Stt) for 0 < (/Jg < tt 
Nx{ps) ± Nx {ps + tt) for n < pa < 2tt . 


The above equations provide sums, (pa) and 
N^{pa), and differences, Npp{pa) and Np^{pa), of counts 
depicted in Fig. [T] (b). While the sums are constant, the 
differences oscillate around zero, and the asymmetry. 


e{ps) = 


N^{pa) + N^iPs) 


3 ctodA® - o-ouA® . _ 

= xPi - —^ sm(^g + p) , (7) 

^ O-QD + O-QU 

in the range (/3g S [0, 27r) has the functional form 


e(v3s) = esin((/jg-b (^), (8) 


independent of beam intensity and target density. Since 
the spin coherence time (SCT) of the in-plane vector po¬ 
larization is long (tsct ~ 300 s), the polarization is 
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FIG. 2. Measured asymmetry e(</3s) of Eq. (O fitted with 
e{ipa) of Eq. ([Sll to extract amplitude e and phase ip, using 
the yields of Fig. [T] (b) for a single turn interval of 

An = 10® turns at a measurement time of 2.6 s <t< 3.9 s. 

assumed to be constant over the duration of the turn 
interval An (1.3s). 

In every turn interval, the parameters e and ip of 
Eq. ([8]) are fitted to the measured asymmetry of Eq. 0. 
An example is shown in Fig.[^ The procedure is repeated 
for several values o f in a certain range around Vg = 'jG 
(see e.g., Fig. 5 of [l3|l. 

A fixed common spin tune Inf’ll = 0.160975407 is cho¬ 
sen such that the phase variation <p>(n) is minimized, as 
shown in Fig. [3] (a). The spin tune as a function of turn 
number is given by 

|j..(n)| = Izzf I + + Az.s(n), (9) 

independent of the particular choice of because a 
different choice for is compensated for by a corre¬ 
sponding change in Ai'g{n). 

Without any assumption about the functional form of 
the phase dependence in Fig. [3] (a), one can calculate 
the spin tune deviation Ai's{n) from by evaluating 
dip{n)/dn using two consecutive phase measurements, 
corresponding to a measurement time of 2.6 s. At early 
times {aip ss 0.06, see Fig. [3] (a)) the statistical accuracy 
of the spin tune reaches = 1.3 x 10“®, and towards 
the end of the cycle « 0.15) = 3 x 10“®, due to 

the decreasing event rate. 

An even higher precision of the spin tune is obtained by 
exploiting the observed parabolic phase dependence, fit¬ 
ted to (p{n) in Fig. [3] (a), which indicates that the actual 
spin tune changes linearly as a function of turn number. 
As displayed in Fig.[3](b), in a single 100 s long measure¬ 
ment, the highest precision is reached at t « 38 s with an 
error of the interpolated spin tune of = 9.7 x 10“^^. 

The achieved precision of the spin tune measurements 
agrees well with the statistical expectation. The er¬ 
ror of a frequency measurement is approximately given 
by (7/ = \/&/Nj( ttsT), where N is the total number 


time t [s] 

0 20 40 60 80 



FIG. 3. (a): Phase as a function of turn number n for all 
72 turn intervals of a single measurement cycle for = 

0.160975407, together with a parabolic fit. (b): Deviation 
Aus of the spin tune from as a function of turn number in 
the cycle. At t « 38 s, the interpolated spin tune amounts to 
\vs\ = (16097540628.3 ± 9.7) x 10"“. The error band shows 
the statistical error obtained from the parabolic fit, shown in 
panel (a). 


of recorded events, i « 0.27 is the oscillation ampli¬ 
tude of Eq. ([S]), and T the measurement duration. In 
a 2.6 s time interval with an initial detector rate of 
5000 s“^, one would expect an error of the spin tune of 
= cr/^//rev ~ 1 X 10“®, and, during a 100s mea¬ 
surement with N Ki 200000 recorded events, an error of 

CT,, « 10-10. 

The new method can be used to monitor the stability of 
the spin tune in the accelerator for long periods of time. 
As shown in Fig. 01 the spin tune variations from cycle 
to cycle are of the same order (10“® to lO”®) as those 
within a cycle [Fig. [3] (b)], illustrating that the spin tune 
determination provides a new precision tool for the inves¬ 
tigation of systematic effects in a machine. It is remark¬ 
able that COSY is stable to such a precision, because it 
was not designed to provide stability below Ri 10“® with 
respect to, e.g., magnetic fields, closed-orbit corrections 
and power supplies. Presently investigations are under¬ 
way to locate the origins of the observed variations in 
order to develop feedback systems and other means to 
minimize them further. 

Several systematic effects that may affect the spin tune 
measurement are briefly discussed below. Terms with a 
vertical vector and a tensor polarization have been omit¬ 
ted in the derivation of e{ips) [Eq. ([7])[. A detailed analy¬ 
sis taking these terms into account shows that py has no 
influence on the spin tune at all, because the particle en- 
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time t [s] 


FIG. 4. Walk of the spin tune during eight consecutive cycles 
with alternating initial vector polarization (black) and 
P(~ (gray). The third cycle is depicted in Fig. |3] (b) as well. 
Cycles with unpolarized beam that followed the p^~ state are 
not shown. 


semble precesses about the j/-axis; py thus merely dilutes 
the asymmetry e{(ps). Although a small tensor polar¬ 
ization of up to ±0.02 leads to higher harmonics in the 
oscillation pattern from which the spin tune is derived 
(Fig. [5]), these contributions alter neither the location of 
the zero crossings nor that of the extrema, and thus have 
no influence on the extracted spin tune. In addition, a tilt 
of the invariant spin axis or misalignment of the detector 
leads to a modification of the magnitude of the measured 
asymmetry e, but neither effect alters the measurement 
of the precession frequency. 

Effects of time-dependent variations of the in-plane 
polarization, acceptance, and flux were studied using a 
Monte Carlo simulation, for which detector rates were 
generated using Eq. ([5]). The analysis, carried out as¬ 
suming these quantities to be constant, showed that even 
extreme variations such as a complete loss of polarization 
or acceptance during a 100 s measurement, does not af¬ 
fect the spin tune determination down to a level of 10“^^. 

The work presented here can be compared to the 
measurement of the muon precession frequency IUmbmI, 
which was determined in the muon {g — 2) experiment 
with a relative precision of Ri 10“® per year [30|. This 
corresponds to an absolute precision of the spin tune of 
Ri 3 X 10“® per year. The higher precision achieved 
here is mainly attributed to the much longer measure¬ 
ment time of 100 s compared to the measurement time of 
600 fis in the muon [g — 2) experiment. Ring imperfec¬ 
tions introducing MDM rotations about non-vertical axes 
make it impossible at this stage to use the new technique 
to directly determine the gyromagnetic anomaly G with 
high precision from the measured spin tune. 

Futnre charged particle EDM searches with an antici¬ 
pated precision of 10“^® e-cm can be carried out in frozen- 
spin mode (9l.[lo|. These investigations, however, demand 
a new class of storage rings. Using an existing machine. 


one could perform a first direct measurement of the pro¬ 
ton or deuteron EDM using an rf Wien filter (sil - f^ . In 
this case one has to cope with the fast spin precession due 
to the deflection and focusing in the magnetic elements. 
The precision determination of the spin tune allows one 
to lock the phase of the spin precession to the rf phase of 
the Wien filter, using a feedback system. 

The method to determine the spin tune, described in 
this paper, provides a precision tool to map out field im¬ 
perfections, orbit corrections, and beam instabilities. It 
can be readily extended to protons. In addition, increas¬ 
ing the measurement period by a factor of ten with tqct 
of a few hundred seconds further increases the precision 
of Vs by about the same factor. 

This paper presents the most precise measurement to 
date of the spin tune in a storage ring. The current pre¬ 
cision reaches a level of = 10“^° for a 100 s mea¬ 
surement. The new method will have a huge impact on 
future precision measurements in storage rings, such as 
the determination of electric dipole moments of charged 
particles. 
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